On Metal— Insulator Transitions due to Self— Doping 
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We investigate the influence of an unoccupied band on the transport properties of a strongly 
correlated electron system. For that purpose, additional orbitals are coupled to a Hubbard model via 
hybridization. The filling is one electron per site. Depending on the position of the additional band, 
both, a metal-to-insulator and an insulator-to-metal transition occur with increasing hybridization. 
The latter transition from a Mott insulator into a metal via "self-doping" was recently proposed 
to explain the low carrier concentration in Yb4As3. We suggest a restrictive parameter regime 
for this transition making use of exact results in various limits. The predicted absence of the 
self-doping transition for nested Fermi surfaces is confirmed by means of an unrestricted Hartree- 
Fock approximation and an exact diagonalization study in one dimension. In the general case 
metal-insulator phase diagrams are obtained within the slave-boson mean-field and the alloy-analog 
approximation . 
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I. INTRODUCTION 

One of the most remarkable examples of the failure of 
the independent-electron approximation is the existence 
of so-called Mott-Hubbard insulators. The early transi- 
tion oxides like LaVOg, LaTiOg, V2O3, TiaOg or CrzOg 
show an insulating behaviour, although they contain par- 
tially filled d bands. The reason is the strong Coulomb 
repulsion of the d electrons, which hinders them from 
moving freely. A simple model to describe the situation 
is the Hubbard Hamiltonian with one electron per siteu: 



H 



(1) 



Electrons can hop from one site i to another j via the 
hopping-matrix element . Whenever two electrons (of 
opposite spins) occupy the same site, they repel each 
other with an interaction energy U. 

To obtain insulating behaviour (implying a vanishing 
dc conductivity at zero temperature) an energy gap for 
charge excitations is required. The existence of a con- 
ductivity gap can be inferred from the spectral density 
of states p{oj). Obviously, for large Coulomb repulsion 
and half-filling there is a gap of order U. The transi- 
tion from a metallic {U — 0) to an insulatina, {U large) 
state is still subject of intensive investigationsH. For gen- 
eral lattices (without special symmetry properties) there 
is reason to believe thaX Xhe transition takes place at a 
finite value U = Uc ^ (MM. 

The Hubbard model deals with one orbital per site and 
ignores any additional bands which might be present. As 



was first pointed out by Zaanen, Sawatzky and AUenS 
this is not sufficient for a classification of the transition 
metal oxides. The following Hamiltonian represents a 
natural extension of the Hubbard model by including 
one additional spin-orbital with creation operators c^^ 
for each site i: 

+ E ^-s- + ^ E (/-^- + ^•^•) • (2) 

In the following we will refer sometimes to the additional 
orbital as a ligand orbital implying that it may belong 
to a ligand atom. Note, that the on-site energy of the 
ligand orbitals has been set equal to zero. Therefore the 
ligand orbitals are higher in energy than the / orbitals 
(A > 0). The system is assumed to be quarter-filled, i.e., 
with one electroa (or hole) per site. 

Zaanen et alu employed a single-impurity approach 
to derive a phase diagram for the transition-metal com- 
pounds. Thereby the transition-metal atoms are taken 
to be independent. This implies a neglect of the hopping 
term in (|^). Moreover, indirect hopping processes be- 
tween the / orbitals via the c orbitals are excluded. The 
calculated band gaps in the density of states then approx- 
imately scale with the hybridization V . Beside the Mott- 
Hubbard regime, where the gap is proportional to U , a 
second region with a gap proportional to A was identi- 
fied. For this phase the term charge-transfer insulator is 
used. As a consequence of the single-impurity approach 
a metallic state was predicted when the transition-metal 
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level lies within the free electron band of the ligand host. 
In contrast, a hybridization gap may be obtained when 
taking into account coherence effects of the transitionpj 
metal atoms. Nimkar et aln and Sarma and BarmanEI 
called the insulating phase due to a hybridization gap a 
covalent insulator and derived the corresponding region 
of the phase diagram. The special feature of a cova- 
lent insulator is, that it undergoes an insulator-to-metal 
transition with decreasing hybridization. The gap does 
not scale with V. LnNiOa (with Ln a rare earth ion oth 
than La) and NiS are believed to belong to this clas; 
They can be driven into a metallic state without a sym- 
metry lowering structural transition by varying pressure 
or temperature. 

If we do not neglect the hopping term Uj in (||) , we can 
also imagine the opposite taking place, i.e., an insulator- 
to-metal transition with increasing hybridization. Elec- 
trons hopping into ligand levels create holes in the Hub- 
bard system, which may move and contribute to the 
conductivity. A mechanism of this kind was recently 
proposed for the rare-earth compound Yb4As3cl. Such 
an insulator-to-metal transition due to "self-doping" 
would be an interesting new possibility for forming metals 
with low carrier concentration. Moreover, like the Mott- 
Hubbard insulator, it cannot be understood within an 
independent electron picture. 

To investigate in more detail the influence of an addi- 
tional band on the transport properties of a strongly cor- 
related electron system we suggest the following model: 



nLn! 



1,(7 "ij^c i 

+ vY^(flc,, + \i.c). (3) 



Again, the system is assumed to be quarter-filled, i.e., 
it has one electron per site. The Coulomb repulsion is 
chosen to be greater than the critical value Uc of the 
Hubbard system. Therefore we will call the model a 
self-doped Mott system (SMS)in the following. In com- 
parison to the Hamiltonian (^ we have neglected the 
hopping iij between the ligand orbitals. For realistic sit- 
uations this approximation is not justified, neither for 
rare-earth nor for transition-metal compounds. Note, 
that in systems like Yb4As3 the widths of the correlated 
and uncorrelated bands are equal {tij = Uj ) . The Hamil- 
tonian (^ should be regarded as a basic model contain- 
ing the necessary terms for the self-doping mechanism. 
Moreover, whenever we have included the hopping term 
t 

come back to this point later 

In this paper we report first results of an investigation 
of the system described by (||). In Section || we illus- 
trate in a transparent way some model properties. We 
discuss exact solutions in the limit of zero bandwidth and 
of infinite hybridization. Concerning the metal-insulator 
transition in the Hubbard model two different situations 
have to be distinguished. In the case of a nested Fermi 



jj we found no change in the general behaviour. We will 



surface the Hubbard model at half filling is an insulator 
for any U > (M. We investigate this case for the SMS by- 
means of an unrestricted Hartree-Fock approximationll3 
and present the results in Section IH. Additionally the 
one-dimensional case is studied by exactly diagonaliz- 
ing a finite system. Two early approaches have domi- 
nated the discussion about the Mott-Hubbard transition 
for general lattices without nesting: (i) the Brinkman- 
Rice approachllil, which can be formally derived by using 
slave boson methodsEj, and (ii) the alloy analogy, which 
is contained in the first part of the well-known Hubbard 
III paperllj. We apply these approximation schemes to 
our model and present the results in Section IV . Finally, 
Section ^ contains a short summary and outlook. 



II. THE MODEL 



In the case of vanishing hybridization the SMS reduces 
to a Hubbard model with a separated ligand level. In the 
first column of Fig. |^ we sketch the corresponding den- 
sity of states. For U > Uc three different regions exist 
(indicated by (2), (3) and (4) in that figure): 




(2) 




/ \ 




(I) 



(II) 




(4) 



Ji 



FIG. 1. Schematic plot of the density of state p{u}) for 
the SMS in the extreme limits of vanishing and infinite hy- 
bridization. The occupied (electron-removal) part is hatched. 
With IMT and MIT we indicate an insulator-to-metal and a 
metal-to-insulator transition, respectively. Different param- 
eter regions are shown (see text). 

• Mott-Hubbard insulator: U < A. When adding 
one electron to the system the lowest energy states 
contain an additional double occupied / orbital. 
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The gap is of order U (Fig. |, (2)). 

• Charge transfer insulator: U > A. An additional 
electron avoids the Coulomb repulsion felt in the / 
subsystem and prefers to occupy the ligand orbital. 
The energy scale of the gap is set by A (Fig. |l|, (3)). 

• Due to the hopping term tij the lower and upper 
Hubbard bands possess a finite bandwidth W{U). 
The ground state of the SMS contains occupied 
ligand orbitals if the gain of kinetic energy due 
to holes in the Hubbard subsystem overcomes the 
transfer energy: A < W. Adding or removing one 
electron changes the energy only infinitesimally. 
The system is metallic (Fig. |l|, (4)). 

We turn now to the case of infinite hybridization. For 
F A it is convenient to use a basis in which the hy- 
bridization term is diagonal. Introducing: 



ifia + Cicr) ; &icr — (/it 



we obtain for the Hamiltonian 



H — H„ + Hh + Hn 



C^a) , (4) 



(5) 



with: 



H. 



ah 



A 

T 

u 



i 

- (4 + ««T + ^.t) ( «a + ^l) ( «a + ^a) }• 

For V = oo the a orbitals remain empty and we can 
restrict ourselves to the part Hi,. The SMS reduces to 
a Hubbard model with a hopping matrix tij/ 2 and a 
Coulomb repulsion C7/4. The result is easily understood. 
When the two orbitals at a given site hybridize strongly 
with each other, the eigenstates are given by the bonding 
h and antibonding a linear combinations. Only the h or- 
bitals are populated to which the original / orbitals con- 
tribute by one half. The Coulomb repulsion felt by two 
electrons occupying one h orbital is reduced to a quarter 



and the hopping between two h orbitals is reduced to one 
half. The trivial but important observation is, that the 
Coulomb repulsion is stronger reduced than the hopping. 
Since only the relation between these two quantities is 
important, the SMS describes for U/2 < Uc a metal and 
for U/2 > Uc an insulator (V = oo). This is indicated in 
the second column of Fig. |l| by (I) and (II) , respectively. 

Having discussed the extreme limits, two transitions 
seem to be possible when changing the hybridization V 
from zero to infinity: For A < W{U), U > 2Uc a metal- 
to-insulator and for A > W{U), Uc < U < 2Uc an 
insulator-to-metal transition should take place with in- 
creasing V . Here W{U) is the effective band width of the 
lower Hubbard band. The two possibilities correspond to 
the aforementioned covalent insulator and a self-doped 
metallic state, respectively. It is worthwhile to note that 
the second scenario is realized only if Uc > 0, i.e., for a 
finite value of the critical Coulomb repulsion Uc in the 
corresponding Hubbard model. 

In order to gain better insight into the behaviour of 
the system for finite hybridization we have considered 
the zero bandwidth limit of the SMS: 



H 



/ 



U n^n^^ 



h.c. 



(6) 



The Hamiltonian is easily diagonalized, yielding the 
eigenvalues Em ^ and eigenstates |to'"-' ) for n electrons 
(n = 1, 2, 3). The / density of states is given by: 



'^M= E{ 



X5(c. + i?r)-i?(f= + l)) 



c<5(.. + £;(f=-i)-£;r))}. 



(7) 



Here Nc = 1 and the sum is only over one site. Fur- 
thermore, I 0'-^'^) denotes the ground state. A similar 
expression holds for the density of states of the c elec- 
trons. The total one is given by p{uj) = p-^{uj) + p'^{uj)- 

In Fig. H we show the results of the diagonalization. 
Plotted are the positions of the poles and their weights 
in the density of states. For small hybridization we can 
clearly identify the lower and upper Hubbard peaks, as 
well as the peaks caused by adding one electron to the lig- 
and orbital situated at « (/ f^c^) and w « A 



f^c'^). Due to the different effect of the hy- 



bridization on a singlet and a triplet state the / ^c^-peak 
splits up. As a result weight is pushed towards the lower 
Hubbard peak. We stress that there is no direct weight 
transfer between the lower and upper Hubbard peak. In- 
stead, the weight of the lower Hubbard peak is always 
equal to unity. Thus a simple self-doping picture which 
predicts a metallic state for an infinitesimal small hy- 
bridization V does not hold true in that case. Note, that 
the gap E^gap = 



(2) 



i^g"' — 2E'^''' between the occupied 
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and unoccupied part of the spectra continously decreases 
from iJgap = A to iJgap « t//4. When the fuU Hamilto- 
nian is used the peaks will mark the position of centers 
of bands. An overlap of the bands accompanied by an 
insulator-to-metal transition is expected, when the hy- 
bridization V exceeds a critical value Vr. 



FIG. 2. Schematic representation of the excitations from 
the ground state of the SMS in the zero bandwidth limit on 
adding and removing an electron. Indicated are the weights 
in the spectral density. Various values of the hybridization 
are used (A = 2, [/ = 4). 



III. NESTED FERMI SURFACE 

We turn now to a treatment of the full Hamiltonian (|^) . 
In this section we consider lattice symmetries and hop- 
ping matrices which lead to a nested Fermi surface. The 
critical Coulomb repulsion in the corresponding Hub- 
bard model vanishes. In one spatial dimension this is 
the general behaviour. The dimension d = I allows for 
an exact diagonalization study. Higher dimensions are 
treated within the unrestricted Hartree-Fock approxima- 
tion, which getting exact in the limit [/ — > and in the 
atomic limit U = oo. Note, that in the case oi Uc = 
small Coulomb repulsions are the most interesting. Af- 
ter the discussion of Section || the behaviour for values 
U ^ Uc should be representative also for all U > Uc- 



Here N is the number of sites in the lattice. It is conve- 
nient to introduce the density of states of free particles: 



Po(w) ^-J^^^i^ 



(9) 



In the calculations presented in this paper it suffices to 
specify Po(j-^) instead of the lattice structure and the hop- 
ping matrix elements. In order to proceed analytically as 
far as possible, we will use in the following a rectangular 
density of states: 



Po(w) 



1 

2W 



Q{W -uj)e{W + uj), 



(10) 



where Q{uj) is the step function. This choice maintains 
the antiferromagnetic instability for arbitrarily small val- 
ues of U and represents the situation of a nested Fermi 
surface. 

In the mean-field approach the interaction term of the 
Hamiltonian (H) is linearized, i.e., 



/ / 



{nl) + {n{^)nl-{n{,){nl). (11) 



'il ■ 



The electrons interact via a molecular field. Depending 
on the anticipated symmetry of the ground state, differ- 
ent ansatzs for the expectation values ( n£. ) are made. 
In the case of antiferromagnetic order two sublattices A 
and B are introduced. The following ground states are 
considered 



• paramagnetic 

• ferromagnetic 

• antiferromagnetic 



— I 



A) 



i(l -(5 + m), 
\{l-5 



m) 



m] 



n; 



n 
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These are the main states in the phase |diagram of the 
doped Hubbard model calculated by PeniM. We indicate 
now the construction oLthe phase diagram for the SMS. 
Following Langer et alrB we solve the equation of motion 
for the one-particle Green's function ({f^^. ; /J ))^ and 

((cg^; ci ))^. From the imaginary part the spectral 



densities pi (oj) and pi (uj) are obtained. For zero tem- 
per ka 

perature the magnetization m and the self-doping 5 are 
given by the solution of: 



A. Unrestrictive Hartree-Fock approximation 

We assume in the following a bipartite lattice. Ad- 
ditional lattice symmetry properties enter the calcula- 
tion only via the band dispersion e-^, which is the Fourier 
transform of the hopping matrix 



1 



k 



(8) 



-y 

k,a — c 



duj = 1 



/ pi^i^)doJ = lil-S + m) (12) 

k — C30 
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where we have introduced the chemical potential fi. The 
equations (^2|) are solved numerically by iteration using 
for po(<^) a rectangular form. Knowing the parameters 
TO, S and fj, for the ground state of a given symmetry the 
energies of the various states are calculated. The state 
with the lowest energy is the true ground state and de- 
termines the spectrum. 

First we discuss the result for V — 0. In the antifer- 
romagnetic phase the model corresponds to an insulator. 
The corresponding spectral density is the one indicated 
in Fig. |l|, (2) and (3). For A > Vl^ the ground state is an- 
tiferromagnetic for all U. This phase can also be present 
if A < W. In particular we obtain an insulating antifer- 
romagnetic solution for small t/ < A up to A = 0. This 
reflects the fact that the ground state does not contain 
occupied c orbitals. The / subsystem is undoped. 

We now turn to the results for finite hybridization. 
Our calculations show that the insulating phases survive, 
when the hybridization V is changed from zero to a fi- 
nite value. In the limit V ^ oo the ground state is an 
antiferromagnetic insulator for all values of A and U. No 
insulator-to-metal transition takes place. In Section || 
we have argued that this behaviour requires a vanish- 
ing critical Coulomb repulsion for the Hubbard model. 
Therefore the obtained result confirms our prediction. 

As mentioned before, in the limit V oo the system 
is always insulating. This finding implies a metal-to- 
insulator transition due to the opening of a hybridization 
gap in a certain parameter regime. We have identified a 
regime with Vc = and one with a finite critical value 
of the hybridization for the transition to an insulator. 
Details will be presented elsewhere. Here we character- 
ize the two regimes only roughly by A < J7 < 2W and 
2W <U < 4W, respectively (A < PF). 



B. Exact diagonalization in one dimension 

In one spatial dimension the critical Coulomb repulsion 
Uc of the Hubbard model with nearest neighbour hopping 
vanishes. The model should be insulating for all values 
of U and V { A > W). Extending the Hamiltonian (||) 
to (H), i.e., including a hopping term between different 
c orbitals, does not change the result qualitatively. In 
the limit V ^ oo the Hamiltonian (Q) reduces to a single 
band Hubbard model with a Coulomb repulsion U /4 and 
hopping matrix elements {tij -\-tij)/2. A self-doping tran- 
sition can only be expected within the parameter regime 
A > W{U) +W SiudUc <U <2{1 + W/W) Uc. With 
W and W we denote half of the bandwidth of the / and 
the c bands, respectively. Here the conditions on U are 
less restrictive as before. 

To verify the prediction we have performed an exact di- 
agonalization of the Hamiltonian (|^) on a four site chain. 
Treating the extended model in this approach is no addi- 
tional difficulty. Before presenting the results we first in- 
dicate the calculations. The one particle excitation spec- 



trum is defined by Eq. (|7|). The first and second term give 
the electron addition and electron removal spectrum, re- 
spectively. The spectral density p{lu) for small clusters 
can be obtained bjt,the Lanczos algorithm in form of a 
continued-fractionll3. In order to obtain the electron re- 
moval spectrum, we start from the vector ipQ — fi^ V'o^'^' 
and then tridiagonalize the Hamiltonian by the Lanczos 
method. The positions of the poles and their weights 
follow from diagonalizing the tridiagonal matrix. The 
electron addition spectrum is obtained similarly. 
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FIG. 3. Evolution of the spectral density of the extended 
model ^ in one dimension. Results are shown for a four 
site chain. Periodic and antiperiodic boundary conditions are 
used. The parameters are given by: U = 4, A = 2, W = 0.25, 
W = 1.5. The solid and dashed lines indicate the / and c spec- 
tral density, respectively. The vertical line marks the position 
of the chemical potential. The peaks have been broadened by 
e = 0.05. 

Results are shown in Fig. |3| We follow the 
usual representation and plot the (S-function of 
equation (|^) with a (arbitrary) width e — 0.05 

{S{^-Up,i) ^ A pseudo-gap 

seems to exist between the electron removal and the elec- 
tron addition spectra. This is caused by the broaden- 
ing. In all our calculations there is clearly a gap between 
the occupied and unoccupied part of the spectrum. The 
energy regime without weight is much larger than the 
typical distance between the poles. To summarize, the 
Lanczos calculations predict an insulator for the extended 
model (A > VF -I- W) in one spatial dimension. Thus 
in the extended model the self-doping transition is re- 
stricted to a limited parameter regime. This behaviour 
is similar to the one of the SMS. 
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IV. GENERAL LATTICES 

In this section we treat the Hamiltonian within (i) 
the slave-boson mean-field and (ii) the alloy-analog ap- 
proximation. We restrict ourselves to a paramagnetic 
ground state. Nevertheless correlation effects are in- 
cluded and both approximations predict a finite criti- 
cal Coulomb repulsion for the Mott-Hubbard transition. 
The restriction to a paramagnetic ground state does not 
allow the resolution of details of the underlying lattice 
symmetry. The results are therefore interpreted as the 
behaviour of the SMS on a general lattice. We have ex- 
cluded here the case of a nested Fermi surface. Again 
it suffices to specify Po{^)- In the slave-boson approach 
we will use a rectangular form as before and in the alloy 
analogy we will choose an elliptical one: 



(13) 



A. Slave— boson mean— field approximation 

In the slave-boson approach four bosonic fields are in- 
troduced with creation operators d/, s}^ and e} . They 
correspond to the four possible states of an / orbital, i.e., 
double occupation, single occupation with spin cr = ±1 
and empty state. The following constrains ensure that 
only the physically relevant part of the enlarged Hilbert 
space is considered 



^la-^ta ^ "i Jia-Jia 



(14) 



Having introduced the bosonic fields the Hamiltonian (|^) 
projected onto the physical subspace can be expressed in 
the form 

i 

(15) 



+ F^(z;„/,U.+h.c, 



where z^^ = e^'Sj^. + s-_^d^. As long as the constrains 



(O) are satisfied, the following replacement is possible: 
1 



1 - d/d, - sls^^ 



1 ^i-a^i-a 



(16) 



Equations (|1J), (|15|) and ( |16| ) provide an exact repre- 
sentation of the projected original Hamiltonian. In the 
mean-field approximation the different representations 
Zj^^ and are not identical any more. The form chosen 



in (|T^) gives the correct result for the Hubbard model in 
the fimit U = 0. 

In order to perform the mean-field approximation in 
the paramagnetic phase, we have replaced the operators 
by real, spin- and site-independent expectation values e, 
s and d. Due to the constrains only one of these values 
is independent, say {= b). The Hamiltonian can then 
be written in the form: 

= - A E + E + 

+ yz^(/,tc,.+h.c.), (17) 



where 



inf-b){Vb+^l + b-2nfy 



— / \ — 



(1 - nf)nf 

i| ) ) . The mean-field approach yields 
an insulating state for vanishing z = 0. In that case the 
renormalized hopping and hybridization vanish. All elec- 
trons occupy / orbitals and are localized. To calculate 
the unknown parameters Uf and b we proceed in a simi- 
lar way as described in Section ( HI A ). The equations of 
motion are solved for the one-particle Green's functions 
((/c ; /- » ,((cr ; c! )) and ((cr: ; /J )) and the 
corresponding spectral densities are calculated. The fol- 
lowing conditions must be obeyed: 



^E / {pi>)+pi»)d-=i 

^E / Pli-)du;^2nf (18) 



k,<T — OO 



k,(T — OO 



dH . 
'db 



db 



E / Hp{Md^ + u^ 



k,(7 —OO 



k.cr — c 

We have solved equations ( p^ ) by using for po{uj) a rect- 
angular form. An analytic expression for the surface 
z = in the {U,A,V) space is obtained, starting from 
the metallic regime. 

For = we— recover the well known result of 
Brinkman and Ricdlil. The critical Coulomb repulsion for 
the Mott-Hubbard transition is found to be Ubr = . 
The value is independent of A even for A < W , because 
the renormalized bandwidth reduces to zero in the insu- 
lating phase. For U < [/br a metallic state exist for all 
values of V. For U > Ubr an insulator-to-metal tran- 
sition is obtained at a critical value of the hybridization 
given by: 



6 



( ^4A2 + U{U - m) - 2A' 



(19) 



The resulting phase diagram is shown in Fig. ^ To sum- 
marize, the slave-boson approach supports the existence 
of an insulator-to-metal transition with increasing hy- 
bridization. In contrast to the arguments given in Section 
H the transition is also found for U > 2Ubr- The transi- 
tion takes place at a finite critical value of the hybridiza- 
tion which in the limit [/ — )■ oo is given by Vc = A/\/2. 



20.0 



12.0 - 




FIG. 4. Phase diagram for the SMS in the slave-boson 
mean-field approximation. The solid lines separate insu- 
lating and metallic regimes. The dashed line indicates the 
Brinkman-Rice value for the critical Coulomb repulsion. Be- 
low the line the system is metallic for all values of A and V 
{W =1). In the limit U ^ oo a metallic phase exists for 
V>Vc^ A/^/2 . 

We have also applied the slave-boson mean-field ap- 
proximation to the Hamiltonian (||). For A > 1^ the 
result is similar to the one obtained for the SMS. As be- 
fore we denote with W half of the bandwidth of the c 
band. The expression for the critical hybridization be- 
comes here: 







Vr 



a ^ia^P'^+U{U-iW)-2al3 
U 2 



otherwise 



where a = A (A — W) 



W 



(20) 



As before in 



the SMS, a self-doping transition takes place for all values 
U > Ubr. 



B. Alloy— analog approximation 

As first pointed out by Velicky et a/Jli the "scatter- 
ing correction" , described ir£3, can be viewed as an alloy 
problem. Applying tiifi approximation to the SMS we fol- 
low the work of Ref.E£l for the periodic Anderson model. 
Alternatively, a Green's function decoupling scheme in 
the spirit of the "scattering correction" can be applied. 
Both approaches are equivalent. The alloy-analog ap- 
proach starts from an intuitive physical picture. An elec- 
tron with spin a can only hop onto the / orbital situated 
at site i, if the orbital is either empty or occupied by 
an electron with spin —a. The electron with spin a can 
therefore be thought of as moving in a static random 
potential with eigenvalues and probabilities given by 



—A , with probability I — 
-A + U , with probability {n{_ 



(21) 

The many-body Hamiltonian is replaced by a one 
particle one with disorder and is of the form: 

+ ( flc,, + h.c.) . (22) 



In the following we assume site- 
expectation values in ( pT[ ) {uf 



and spin-independent 



'it ' 



The Green's function G corresponding to the Hamilto- 
nian (|2^ ) has to be averaged over all possible configura- 
tions of the random potential which can be considered 
being due to impurities. The averaging cannot be per- 
formed exactly. To solve the alloy problem the coherent 
potential approximation (CPA) is used. The averaged 
Green's function G follows from an effective Hamiltonian 
containing a dynamical mean-field S(ix>): 

E(cj) will be determined in (p6|). For an elliptic density of 
states po{uj) the averaged / Green's function is therefore 
of the form: 



W2 



UJ 




(24) 



The other Green's functions, G-^'^ and G'^'^, can be ex- 
pressed in terms of G^-^ . A scattering matrix T is intro- 
duced for each configuration via: 
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G = G + GTG. (25) 

The dynamical mean-field is determined by demanding 
that the scattering matrix vanishes on average: T = 0. 
This yields an expression for Tj{uj) of the form: 

^ E- (-A - EH) (-A + U- EH) G^-^(w), 

(26) 

where E — —A + U nj. Eliminating E(a;) from ( p^ ) and 
( p6| ) gives a cubic equation for G^^{uj). By guessing a 
starting value for n j one solves this equation and calcu- 
lates the / and the related c density of states p ^ {uj) and 
p-^(w), respectively. An improved value of is obtained 
from the solution of: 

J p^{uj)duj = nf 
^ (27) 

— oo 

The equations are iterated until self-consistency is ob- 
tained. 
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FIG. 5. Evolution of the spectral density with increasing 
hybridization for the SMS in the alloy-analog approximation. 
The filled area indicates the occupied part of the spectra. The 
used parameters are given by: U = 4, A = 2, W = 1. 

It is well known that the alloy analogy gives a critical 
Coulomb repulsion of Uaa — W for the metal-insulator 
transition in the Hubbard model. A typical result for 
U > Uaa is shown in Fig. |[ The total weight of the 
lower Hubbard band is given by 2 (1 — n/), independent 
of the on-site energy of the ligand orbitals. For any fi- 
nite value of the hybridization, Uf is less than 1/2 and the 
system is metallic. This corresponds to the simple self- 
doping argument mentioned in Section | and The ob- 
tained weight transfer from the upper to the lower Hub- 
bard band is not surprising, since the expectation value 



n f determines the probability for the two on-site energies 
—A and —A + U of the / orbitals in the alloy problem 
(Eq. (pT|)). The comparison of Fig. || and || reveals, that 
the overall evolution of the spectrum is well described 
within the alloy-analog approximation. The exact peaks 
obtained in the limit of zero bandwidth are situated in 
the middle of the alloy bands. The only exception is the 
band which is the singlet /^c^-peak in the zero band- 
width limit. It is shifted toward higher energies. In the 
limit W ^ and V ^ oo the gap between the lowest 
two peaks is given in the alloy approach by C//2 and not 
by [//4 as in the exact result. 



V. CONCLUSION 



We have investigated the effect of a hybridization V 
between an orbital with strong Coulomb repulsion and 
a ligand orbital on the transport properties of a Hub- 
bard system. The Coulomb repulsion U was taken to 
be larger than the critical interaction Uc for the Mott- 
Hubbard transition of a system without ligand orbitals, 
i.e., U > Uc- The ligand orbitals lay higher in energy 
than the correlated ones (A > 0) and the occupation is 
one electron per site. 

For zero hybridization two regimes can be distin- 
guished. Depending on the on-site energy of the ligand 
orbitals the ligand level enters the lower Hubbard band 
or not. In the first case a metal-to-insulator and in the 
second case an insulator-to-metal transition takes place 
with increasing hybridization. Both forms may exist in a 
range of the Coulomb interaction U. The former occurs 
for t/ > 2 Uc and the latter for [/< 2 C/c. 

As known from the exact solution of the Hubbard 
model the critical Coulomb repulsion Uc vanishes in one 
dimension. We confirmed the absence of a self-doping 
transition in that case by an exact diagonalization study. 
In higher dimensions lattice symmetries and hopping ma- 
trices which lead to a nested Fermi surface also give 
L^c = in the corresponding Hubbard model. In these 
cases the self-doping transition again should not exist. 
We confirmed the prediction by means of an unrestricted 
Hartree-Fock approximation. Moreover, we shortly dis- 
cussed the opening of a hybridization gap in the regime 
U > 2Uc and A < VF. 

In the case of C/c ^ the resulting picture is less con- 
clusive. Both, the slave-boson and the alloy-analog ap- 
proach confirm the existence of a self-doping transition. 
A restriction to a limited U regime is not found here. 
Moreover, the transition from a metal to an insulator for 
values U > 2Uc and A <^W does not occur. 

One should appreciate that one is dealing here with 
a rather intricate problem. The number of different ap- 
proximation schemes for investigating the Mott-Hubbard 
transition for finite Coulomb repulsion is rather lim- 
ited. To explore the possibility of a self-doping transition 
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one must be able to deal with the difficult parameter 
regime U Uc- Moreover, in addition to the many- 
body Coulomb term also the hybridization term has to 
be treated. The approximations applied here do not con- 
tradict the existence of an insulator-to-metal transition 
with increasing y in a restrictive parameter space for U. 



At present it is not clear which of the findings are due to a 
particular approximation made, i.e., slave-boson mean- 
field or alloy analogy and which ones are independent 
of it. To resolve this problem further investigations are 
required. 
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